Solid-state NMR quantum computer with individual 
access to qubits and some its ensemble developments 



> 
o 

m 

(N 



K. A. Valiev and A. A. Kokin 



Institute of Physics and Technology, Russian Academy of Sciences, Nakhimovskii pr. 
36/1, Moscow 117218, Russia, E-mail: qubit@mail.ivvs.ru 

Abstract 



Here we made an analysis of the principles of a semiconductor NMR quantum 
computer and its developments. The known variant of an individual-access computer 
| (B. Kane) and alternative solid-state bulk-ensemble approach versions allowing to avoid 

some difficulties in implementing the first variant are considered. 

o 

o\ . 

O Introduction 

On 

Atomic nuclei with a spin quantum number 1 = 1/2 seem to be the natural candidates for 
qubits — two-level quantum elements in quantum computers. The early NMR quantum 
computers approach suggested in 1997 by two teams of researchers independently pi and 
then confirmed in experiments || ||. In this approach it was used several diamagnetic 



organic liquids whith individual molecules, having a small number of tied non-equivalent 
nuclear spins-qubits, acting as almost independent quantum computers. Such an effect is 
due to the fact that vigorous rotational and translational Brownian motion of molecules in 
a liquid to a great extent averages both intra- and intermolecular dipole-dipole nuclear spin 
interactions. Decoherence time of the spin states thus turns out to be relatively large (several 
seconds or more). Only scalar intramolecular interaction between nuclear spins remains not 
averaged. It is described by the Hamiltonian X)i<jiylilj where Zy is an interaction constant. 
In particular, for the chloroform molecule 13 C 1 HCl3 this interaction is essential only for 
proton 1 H and carbon isotope 13 C and it is this interaction around which unitary quantum 
operations on qubits have been performed. 

Since in liquids the nuclear spins are weakly coupled with the environment, we may 
restrict our consideration to nuclear spins of an individual molecule (reduced quantum en- 
semble) rather than deal with a huge number of nuclear spins of all molecules of the liquid. 
For finite temperatures, the reduced ensemble is in the mixed state described by a diagonal 
density matrix. The nonzero elements of this matrix are the occupancies of associated spin 
states. For such a system could perform quantum computation, it should be initialized prior 
to data entry. Initialization means the separation of lower- dimension blocks from the density 
matrix of the initial mixed state that must have properties similar to those of pure states (a 
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pure-state diagonal matrix may have only one nonzero element). States described by such 
matrix are called effective or pseudo-pure. All unitary quantum computations are performed 
using these states. Several approaches to initialization have been reported |2], £|, |5|] . 

In an ensemble computer, many molecules-minicomputers, being nearly independent 
one another, act in parallel and can thus be controlled by operations on the entire (macro- 
scopic) volume of the liquid. These operations are well-known in high-resolution NMR tech- 
niques. Access to individual qubits is replaced by simultaneous access to related qubits in 
all molecules of a bulk ensemble. Computers of this type are called bulk-ensemble quantum 
computers. They, in principle, can operate at room temperature. However, for magnetic 
fields used in conventional NMR spectrometers (no more than 20 T) relative spin polar- 
ization, specifying signals from pseudopure states, is very small even at low temperatures. 
In addition, the NMR signal intensity exponentially drops with the number of qubits ||. 
Therefore, an organic-liquid quantum computer operating at room temperature cannot have 
much more than 10 qubits in a molecule. Estimates show, however, that a quantum com- 
puter of practical value, i.e. with the number of degrees of freedom (or the dimension of 
the Hilbert space of a set of qubits) larger than in conventional customary computers, must 
have more than 40 qubits (2 40 ~ 10 12 conventional bits). 

Another radically new and still unimplemented design was proposed in || . It involves the 
formation of an artificial multiple-spin system with access to individual nuclear spins-qubits. 
It is suggested to use an semiconductor MOS structure on a 28 Si spinless substrate into a 
thin layer of which 31 P stable phosphorus isotopes, acting as donors, are implanted. These 
donors substitute for silicon atoms at the lattice sites, producing shallow impurity states 
with a large effective Bohr radius, and have a nuclear spin number I = 1/2. The number of 
donors or qubits in such artificial "molecules" can be very large. 

In this work (see also our article [JJ]) we analyze in details the principles underlying NMR 
quantum computer operation (discussed in ||) and give a some further development. We 
call attention to the advantages and disadvantages of the early version. Possible alternative 
ensemble designs of solid-state-based nuclear-spin quantum computers that exploit the same 
principles but are free of the above disadvantages are next discussed. 

1 Basic requirements for semiconductor structures 

In the discussed Kane's approach temperatures must be low enough for electrons of the 
donors occupy only the lower spin state in a magnetic field; i.e. T <C I^Bjk where 
fiB = 9.27- 10~ 24 J/T is the Bohr magneton, B is the induction of an external magnetic field 
and k = 1.38 ■ 10~ 23 J/K is the Boltzmann constant. For B < 2 T we have T > 0.1 K which is 
much lower than the temperature for freezing-out donor electron states. Hence, the donors 
will remain for long in the neutral ground orbital S state (D°— state). 

Each donor atom having a nuclear spin must be fairly accurately located under "its" metal 
gate (gate A) separated from the semiconductor by a thin insulator film (for example, several- 
nanometer-thick silicon oxide). Gates A form a linear array of arbitrary length and period 
I (one artifical "molecule") (Fig. 1). With the aid of the electric field generated by gates A 
one can affects the electron density distribution near the nuclei in the ground state, thereby 
individually adjusting the resonance frequency of every nuclear spin defined by external 
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magnetic field and hyperfine nuclear-electron interaction. Thus, quantum operations by 
selectively applying resonant radio-frequency (RF) pulses to the nuclear spins of particular 
donors become feasible. 

In deciding on the location of the donor atoms under the gate, one should bear in mind 
that silicon surface is not perfectly smooth; it will always have irregularities like valleys and 
ridges. Accordingly, the near-surface electric field under the gates at depths on the order of 
several lattice constants (in silicon, the lattice constant is 0.54 nm) will also be randomly 
non-uniform. Therefore, the electric field around the donors can be effectively controlled if 
they are more then 10 nm away from the silicon surface. 

The characteristic sizes of the semiconductor structure (Fig. 1) are in a nanometer range. 
Such structures are fabricated with the use of modern nanotechnology techniques, such as 
epitaxial growth, ultra-high-vacuum (UHV) scanning probe nanolithography based on tunnel 
or atomic-force microscopes (AFM) || and electron-beam or X-ray lithography. 

The interaction of the nuclear spins of the donors with the environment can be precluded 
if the silicon and silicon oxide are enough purified of 29 Si isotope. This isotope has a spin 

1 = 1/2 and is contained in amounts of 4.7% in natural silicon. For the number of silicon 
atoms per cm 3 about 5.0-10 22 and the characteristic sizes depicted in Fig. 1, one donor atom 
occupies a volume of 20 3 nm 3 = 8 • 10 _18 cm 3 . On the average, less than one 29 Si isotope 
atom must be present in this volume. Hence, the isotopic purity of silicon should be as high 
as ~ 2 • 10" 4 %. 

In general, III-V semiconductors are inapplicable to quantum computation at this ap- 
proach, since they do not possess spinless isotopes. However, we assume that spinless semi- 
conductor structures can be built around other Group IV elements. Natural germanium has 
only one spin isotope 73 Ge (/ = 9/2) in amounts of 7.76%. Natural carbon has spin isotope 
13 C with a spin number / = 1/2 in quantities of 1.1%. Therefore, spinless structures can be 
made of Ge, Si/Sii_ x Ge x , and SiC materials if they are appropriately purified. 

The nuclear spins of 29 Si and 13 C isotope atoms, which produce shallow impurity states 
in the above structures, can also be viewed as candidates for spins-qubits. 

2 Electron— nuclear spin system of a donor in a mag- 
netic field 

Nuclear and electronic spins of donor atoms in the ground orbital state \l/o( r ) interact by 
hyperfine and dipole-dipole magnetic interactions. After averaging over the orbital electronic 
state, the latter vanishes, and only hyperfine interaction remains, whose Hamiltonian has the 
form (Fermi formula) 

Hjs = MS, (1) 

where 

A=^|vI/ (0)| 2 2 / i B ^ N/ i N -^ (2) 

3 47T 

is hyperfine interaction constant, = Am - 10 _1 T 2 cm 3 /J, /xn = 5.05 ■ 10 -27 J/T is the nuclear 
magneton, and = 2.26 for 31 P 1 ' ) . With the experimentally found constant of hyperfine 
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interaction A = 7.76 • 1(T 26 J (or A/(2nh) = 116 MHz) we obtain for the probability of 
an electron being on the donor nucleus (r = 0)|\E'o(0)| 2 = 0.43 • 10 24 cm -3 . 
The electron-nuclear spin Hamiltonian for a donor atom has the form 

H = 2/i B BS - ^ N /i N BI + AIS, (3) 

four energy levels of which are given by the well-known Breit-Rabi formula. For 1 = 1/2 
(the z— axis is parallel to B) this formula is written as 

A A 



E(F, m F ) = -- - WN Bm F ± - yf 1 + 2m F X + X 2 , (4) 

where X = (2/i B + g^fi N )B/A,F = I ± 1/2 = 1, 0, and m F = M + m = ±1, if F = 1 
or nip = if F = (Here M = ±1/2 and m = ±1/2 are z— projections of electronic and 
nuclear spins accordingly). The energy levels vs. X are shown in Fig. 2. 

For the energy of the ground spin state, F = and mp = (M = — 1/2, m = 1/2); 
hence, we obtain 



E(0, 0) = -A/A - (1/2)^(2^ + SN/i N ) 2 £ 2 + A 2 . (5) 

For the next, excited energy state, F = 1, mp = — 1, with the changed nuclear spin state, 
M = — l/2,m = -1/2, we have 

£(1,-1) = A/4- (2/i B -<? N /i N )£/2. (6) 

Thus, the energy difference between the two lower states of the nuclear spin that interacts 
with an electron whose state remains unchanged is described in simple terms (/xb 3> /in for 
A 2 /(2/i BJ B) 2 « 1): 

A 2 

27ihv A = E(l, -1) - E(0, 0) = A/2 + g Nf i N B + ~ A/2 + g Nf i N B. (7) 

For 31 P donor atoms, the first term in (0) exceeds the second one in fields B < 3.5 T. 
In this case, the nuclear resonant frequency u\ depends largely on the stationary local 
magnetic field Bi oc = AS/ (pn/^n), which affects the nucleus via electron spin polarization 
S. The external RF field b±(t), normal to the constant field B, in NMR acts on nuclear 
spin not directly but through the transverse component of electronic polarization -Bi OC) ± ~ 
AS\\b_i_(t) / (pn/^n-B) • This gain effect was first indicated by Valiev in |L0] . 



Now, putting B = 2T, according to 0, we obtain for the nuclear resonance frequency 
i/ A = 92.6 MHz. 



3 Electronic structure of the ground state of a donor 
and hyperfine interaction constant 

Consider the electronic structure of the ground state of a donor atom. The conduction band 
in silicon has six isoenergetic valleys with minima displaced relative to the Brillouin zone 
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center by vectors kj (j = 1,2, .. .6) (|kj = k ) toward three orthogonal axes of the fourth 
order. These axes are the axes of symmetry of the crystalline structure of silicon. 

Let the direction to neighboring donors (x— axis) coincide with one of the axes of sym- 
metry of the fourth order, for example [100], i.e., be perpendicular to the directions with 
whose four of the six energy-valley ellipsoids are aligned. Then, the effective Bohr radius a 
of electrons for these four symmetric valleys in the indicated direction will be substantially 
larger than in any other direction, since this radius is defined mainly by the transverse ef- 
fective mass m t , which is much smaller than the longitudinal effective mass mi. The value 
of the effective Bohr radius specifies the extent of the wave function of a donor electron and 
hence the characteristic scale of the semiconducting structure along 

Note that for the same reasons, in germanium, unlike silicon, the direction to a neighbor 
should be taken along one of six axes of symmetry of the second order, for example, [110]. 
This axis is perpendicular to four axes of the third order (in the Brillouin zone) along which 
the valley ellipsoids are aligned. 

The orbital wave functions of an electron of a donor atom at the silicon lattice site 
are written in terms of tetragonal symmetry group representation Td and are expressed by 



superpositions like [11 



*(r) = ^a j F J (r)^(k j ,r), (8) 
j=i 

where ^(kj,r) is a Bloch function quite rapidly varying within a distance comparable to 
the lattice constant and corresponding to the six kj minima in the Brillouin zone. It is 
normalized to the unit cell volume Q: 
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k h v)\'dv=l. (9) 



For silicon, it was found JTT] that |^(kj, 0)| 2 = 186 ± 18, so a conduction electron is 
strongly localized at the lattice site. The wave functions i*j(r) describe smooth donor-related 
modulation of the Bloch functions and satisfy the Schrodinger equation with effective masses 

/ h 2 d 2 h 2 ( d 2 d 2 \ a 2 \ 

' +^)-^ L --EAF;(r) = , (10) 



2midz 2 1m,t \dx 2 dy 2 J Aire s eor 



where the Zj— axis is parallel to kj and originates at the donor atom; are the energy 
eigenvalues of the donor electron; and e s = 11.9(Si),eo = 8.85 • 10~ 14 F/cm. 

The ground state of a donor electron in silicon belongs to one-dimensional (non-degenerate) 
representation A\ of the tetragonal symmetry group The wave function of this state can 



be represented as the superposition flD with equal weights Oj = y 1/6 pTTfl : 

v&o(r) = .yi76^ (1S) W^(k j ,r). (11) 



For silicon and germanium, it was performed variational computations with the test 
modulating wave function ]TT] 



if s) (r) = (vra^O- 1/2 exp[-(p 2 M 2 + zf/a^% (12) 
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where p 2 = x 2 + y 2 , a t and b t are variational parameters for the ground state similar to the Is 
state of a hydrogen atom (the ground state under consideration turns into the Is hydrogen 
state if m t = m{). The obtained results somewhat different from the effective Bohr radii 
with the transverse and longitudinal masses: 

at(Si) = 2.50 nm, ^ (Si) = 1.42 nm, 

o t (Ge) = 6.45 nm, a,(Ge) = 2.27nm. (13) 

The results for the ground state energy E^o were the following (the values in parentheses 
stand for experimentally found ionization energies Ed): 

£ d0 (Si) = -0.029 eV(-0.045eV), 

£ d0 (Ge) = -0.009 eV(-0.012eV). (14) 

The states of the donor atom with the same function (r) but having different weight 
coefficients for two-dimensional, Ex, and three-dimensional, T\, tetragonal group represen- 
tations are excited and omitted from consideration. 

The theoretically obtained probability of an electron being found in the vicinity of the 
donor nucleus in silicon, 

|^ (0)| 2 = 6(F/ ls) (0))V(kj,0)| 2 = = 0,042 ■ 10 24 cm- 3 (15) 

TTCti (X\ 

is almost 10 times less than the above defined experimental value, and the ionization energy 
is underestimated by a factor of 1.5. This is explained by the weaker coordinate dependence 
of the electrical potential in equation (10) produced by a donor atom in silicon at small 
distances intermediate between the covalent crystallochemical radius of phosphorus, 0.11 nm 
IP^ I, and the silicon lattice constant, 0.54 nm. In this small region, one can put (.Fj (r)) 2 ~ 

(F/ ls) (0)) 2 ^ (7ra 2 bi)-\ Using the experimental value for |^ (0)| 2 , we find (F/ ls) (0)) 2 = 
3.94- 10 20 cm- 3 . 

Consider now a variation of the ground state energy due to a deviation of the potential 
from the point-charge value |T3[ : 



AE d = ^(F^\o)) 2 -^-F 2 , (16) 
3 47re s e 



where r 2 is the average square of the radius of the region where potential deviation takes 
place. If the deviation of the ground state energy AE d = 0.045 — 0.029 = 0.016 eV, we find 
from (0) p) 1 / 2 = 0.42 nm, quite a reasonable value. 

It is noteworthy that, while computing | x I / o(0) | 2 for germanium with the use of correspon- 
dent variational parameters, one should take into account that the centers of isoenergetic 
ellipsoids are the centers of eight hexagonal Brillouin-zone faces; therefore, the sum over j 
will contain only four terms (or four full ellipsoids), and the numerical value of |^(kj,0)| 2 
will be different. For a 31 P atom, the experimental value of the hyperfine interaction con- 
stant is A = 15 • 10 _4 cm _1 = 45 MHz, the singlet electronic state factor g = 1.56, and 
|^o(0)| 2 = 0,22- 10 24 cm- 3 M. 
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4 Hyperfine interaction constant of a donor atom ver- 
sus electric field 



Consider in more details how the external electric field E induced by the gate potential 
affects the hyperfine interaction constant A. We shall assume that the field is aligned with 
one of the axes of the fourth order. The origin, unlike equation (10), is now on the gate 
surface rather than at a donor atom (Fig. 1). 

With regard to the fact that, at low temperatures, intrinsic semiconductor is essentially 
a dielectric and also considering that the thickness of substrate D ^> /a ~ c ^> d (Fig. 1), 
we shall express the electric field strength at a donor atom via the potential V on gate A, 
assuming that the gate is a circular disk of radius a = I a/2 lying on the surface of a semi- 
infinite dielectric. The expression for the electrical potential at a point with coordinates 



(p, z),p — x + y has the form [14 



(p{p, z) = — arctan , 



2a 2 

(17) 



7T \ p 2 + z 2 - a 2 + \(p 2 + z 2 - a 2 ) 2 + Acfiz 2 } 1 / 2 ' 

hence, for a line passing through the gate center p = near the donor atom, we have 
<p(0, z) = ip(0, c) - E c (z - c) + E' c (0, c)(z - cf/2 + . . . , 



where 



<p(0, c) = — arctan -, 

7T C 



E c = E c (0,c) 
E' c = E' c (0,c) 

are the potential, the electric field strength and his gradient at the donor atom. Equation 
(10) should now be supplemented by the perturbed Hamiltonian = —q(p(0,z). 

In the case of the almost homogeneous electric field in the neighborhood of the donor atom 
(contribution of term with E' c (0, c) is small if ca^/a 2 <l,aj = S7ie eE d0 /q 2 — the effective 
Bohr radius) a shift in the electronic density distribution relative to the position of the donor 
atom nucleus causes the same decrease in A, irrespective of the direction of the electric field 
with respect to z— axis. A correction to the wave function E^ ls \c) modulating the ground 
state of a donor (its coordinates are now p = 0, z = c) is defined by the second-order 
expression in the perturbation theory. Non-diagonal matrix elements of the perturbation 
operator between even F^ 1 (r) states and odd excited states F^ 2p \r) and _F^ 3p ^(r), similar 
to the 2p and 3p hydrogen states, are not zero. The energies of these excited donor states in 
silicon |[l| are E 2p = — 0.0109 eV and E 3p = — 0.0057 eV. Since the excited state functions 
F^ 2p ' ) (c) = _F^ 3p ^(c) = 0, the desired correction, nonzero at z — c, is given only by that part 



d<p{0, c) _ 2V 
dz 

dE(0,c) 



7r a 2 + c 2 ' 



(19) 



AV 



ac 



dz 



(a 2 + c 2 ) 2 



71 



7 



of the second-order expression in the perturbation theory proportional to the ground state 
function R {ls) (c) IIT 



AfToo = -F/'-' ( c)iAg' Ko|^Mi (20) 

where m = 2p, 3p, etc. Since the energy differences between the ground and excited states, 
E 2p — Ed = 0.034 eV, E 3p — = 0.039 eV, . . . , 0.045 eV, are close to each other, we can 
estimate this correction through a second-order change in the ground state energy and then 
through the polarizability \ °f the atom in an electric field by taking one of the factors in the 
denominator equal to some mean value 8E ~ 0.04 eV and carrying it out of the summation 
sign: 

AF^(V) ~ -FWM.ilr'ifcM, 
j 1 ] j {) 26E^ E m -E d 

= -\/{25E).{ X El/2).Fl ls \c). (21) 

The polarizability \ of a donor atom in silicon can be estimated by the formula for a 
hydrogen atom [O using the effective Bohr radius a B = 2 nm, that leads to % — 4vre • 
(9/2) (ag) 3 = 4 • 10~ 32 F • cm 2 . Eventually, we have for correction (20) 

AF, (ls) (V) = -1.55 • 10~ 12 ■ E 2 C ■ F/ ls) (c), (22) 



where electric field E c is given in V/cm. From (|22"D, it follows that the relative correction 
is small when E c ^ 8 ■ 10 5 V/cm. For a field-dependent correction to hyperfine interaction 
constant, one obtains to second-order (~ E 2 ) accuracy 

AA{V) j A ~ -3.1 ■ HT 12 ■ E\. (23) 

Substituting E c = 0.25 ■ 10 6 V, found from (0) for c ~ 2a ~ 10 nm, into (||) yields 
Ao4(y)/A = — 0.19V^ 2 (here V in V) or, for a V^— dependent correction to resonance frequency, 

Av A (V) = AA(V) /A-v A = -17.5 • V 2 MHz. (24) 

This disagrees with results in |J where a frequency tuning parameter a = dAp^iV) / dV =3- 
— 30MHz/V at V =^ is given. In general, the linear part of the Av A vs. V dependence 
may be attributed to different causes. One of this is a non-homogeneity of the electric field 
(if cog/a 2 ~ 1), and other is the built-in electric field in the semiconductor. In the first case 
it is necessary to take into account the last member of the expression (|TB| ) which have the 
first-order perturbation theory nonzero value. The second case can take place when the work 
functions of the gate and substrate, which determine flat-band voltage Vfb> are distinguish. 
In this case, one should put V 2 =>■ (Vfb + ^) 2 m the above expressions. If for example, we 
put Vfb ~ 0.6 V, 

Au A (V) ~ -17.5-(V FB + V) 2 = 

= -6.3 - 21.0 • V - 17.5 • V 2 MHz. (25) 
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Note, that, as V and Vfb approach 1 V, the second-order expressions of the perturbation 
theory become poor approximation and corrections with must be included in fl22|) . 

The relative error in the hyperfme interaction constant 5 A/ A due to technological inac- 
curacy (5p) 2 in placing donor atoms in the z— axis is defined from (|17D and (|23|) by the ratio 
~ 2Sp 2 /a 2 for c/a< 1 and ~ 25p 2 /c 2 for c ~ a. The condition 5 A/ A <C 1 then specifies the 
lower limit for the parameters c and a ~ ^a/2. 



5 Electron— nuclear spin system for two interacting donor 
atoms 

The distance / between donors, as well as the period of the semiconductor structure, is 
taken small enough. In this case, the constant J of effective exchange electron-electron 
interaction J(S a Sb) between two hydrogen-like neighbors a and b due to partial overlap of 
their electronic wave functions in the corresponding direction is the most sensitive to the 
field of the gate J. 

The total spin Hamiltonian for such a "molecule" is 

H = 2/i B B(S a + S b ) + JS a S b + AH = H + AH, (26) 

where 

AH = H & + H h = -# N/ u N 5(/ za + J zb ) + AJ a S a + A h l h S h . (27) 

In the absence of external magnetic field, for a positive value of the exchange interac- 
tion constant J > ground state X S of the "molecule" is singlet. In 0, for well-separated 
hydrogen-like atoms, an asymptotic expression that does not include oscillations due to inter- 
ference of the Bloch functions for various valleys was employed. In terms of the designations 
adopted in our work, it has the form fH^ff 

J{1) «1,6- — q — ■ (l/a t f 2 exp(-2//a t ), (28) 
47ree a t 

where at = 3 nm for silicon (in || , this variational parameter was used slightly lower: at = 
2.5 nm). 

According to variational parameters (|l~3|), the extent of the wave function of a donor- 
atom electron, which is specified by at, in germanium turns out to be 2.6 times greater than 
in silicon. Accordingly, the distance at which the desired overlap of the wave functions of 
neighboring donors is achieved also increases. Thus, in the case of germanium, the scale 
of the semiconductor structure can be considerably extended if the axes of symmetry are 
properly oriented. 

We shall omit for a while the interaction with nuclear spins in (|26|), which here is a 
small perturbation, and after diagonalization of the Hamiltonian we arrive at three triplet 

3 ) This asymptotic expression differs from the known Sugiura's Heitler-London approximation, which 
gives a physically impossible positive value at extremely large 1/ at > 49, 5. 
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(S = l,M — 0, ±1) and one singlet (S = 0, M = 0) energy levels for a dual-spin system 
(M = M a + M h — the projection of total electronic spin) 

E (S, M) = J[S(S + l)/2 - 3/4] + 2fi B BM. (29) 

Varying the potential of gates J, located between gates A, one can control the electronic 
density between neighboring donor atoms and thus the overlap of the wave functions of 
electrons localized on neighboring donors a and b, as well as the constant of their exchange 
interaction J and the constant of indirect scalar interaction between their nuclear spins 1^. 

From Fig. 3, it follows that, for J = 2fi-gB, the singlet and one triplet levels are crossed 
as B grows. In this range of the parameter J, even rather weak effects may essentially 
change the structure of the ground electronic state of the "molecule". The relationship 
J (I) = 2/xb-S defines the necessary interdonor distance, which was estimated in ||. For 
B = 2 T, it is I ~ 10 — 20 nm 3> a*; hence, the upper limit for the size of gates A in such 
fields is /a ~ 10 nm. If germanium is used instead of silicon, rough estimates using the above 
values of the variational parameters yield much greater sizes: I ~ 25 — 50 nm and I a ~ 25 nm. 

Consider now hyperfme splitting of " molecular" levels due to the Hamiltonian of pertur- 
bation AH (|2"TD (J ^> A a , Ab). The main interest is the splitting of electronic levels, crossing 
and especially an anticrossing of states near the level crossing point. 

We shall take into account, following on |Tj|, that commutator [S za +S z b+I za +I z b, H] = 0. 
Since the total electron-nuclear spin projection in the magnetic field direction m a + + 
M a + Mb = m + M is conserved, matrix 16x16 that represents the total Hamiltonian, falls 
into five reduced matrices, corresponding to values m+M = 0, ±1, ±2. The four states, some 
of which are used for measurements of the nuclear spin states, correspond to m + M = — 1. 

We shall represent the perturbation Hamiltonian (|27| ) as the sum of the secular and 
nonsecular parts: AH = AH sec + AH nsec , where 

A4cc = -9n^nB(L + L) + (1/2X4,4 + A b / b z)(£ az + S hz ), (30) 

[A4cc,^o]=0. (31) 

The paired nuclear spins J a> b = 1/2 of the atoms a and b may be in three triplet (J = 
I a + Ih = 1, m = m a, + m b = 0, ±1) and one singlet (I = J a — lb = 0, m = 0) states. 

States of the electronic-nuclear spin system will be designated as \S, M; I, m). They are 
eigenstates of the total Hamiltonian if the nonsecular part of hyperfme interaction AH nsec is 
neglected. 

Let us concentrate here on two cases: 

1) For J < 2fi B B, the lowest triplet state (S — 1, M — —1) is the ground electronic state 
(Fig. 3). First-order energy corrections to ( P9"| ) within the perturbation theory depend on 
diagonal matrix elements of AH: 

AE 1 (S, M; I, m) = (S,M;I,m\AH\S,M;I,m) (32) 

In this approximation, from the ground electronic level Eq(1,— 1) only three sublevels 
determined by the secular part of Hamiltonian are split: 

AEi(l,-l;l,-l) = g N fi N B + (A a + A h )/A, 
A£i(l,-1;0,0) = A£i(l,-l;l,0) = 0, (33) 
A^l,-!;!,!) = -g N fi N B - (A a + A b )/4, 
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The states |1, — 1; 1, 0) and |1, — 1; 0, 0) that have the same value m + M = — 1, remain 
degenerated. In the next approximation determined by the nonsecular part of hyperfine 
interaction, the lowest level E (l,—1) in the electronic triplet for J < 2[i B B is split into four 
sublevels. 

2) For J > 2/zb-B, the electronic singlet 1 , 0) is the ground electronic state. In the 
first approximation, the singlet ground electronic level is split into next three hyperfine 
sub levels^ : 

^(0,0; 1,-1) = g N fi N B, 
A£i(0,0;0,0) = A£i(0,0;l,0) = 0, (34) 
A^O^M) = -g N fi N B, 

Among them one state |0,0;1, — 1) have M + m = —1. If nonsecular part of hyper- 
fine interaction is accounted this state participates in an anticrossing process for the state 
|1,-1;0,0). 

Near the crossing point C for the ground electronic states the nonsecular part of hyperfine 
interaction AH nsec has next four nonzero non-diagonal matrix elements between states with 
the same values m + M: 

(0,0;0,0|A/f nsec |l,-l;l,l> = -(A a + A b )/A, 

(0,0;l,-l|A/J nsec |l,-l;l,0> = (A a - A h )/A, (35) 

(0,0;1,-1|A/J nsec |l,-1;0,0> = (A a + A h )/A, 

(0,0; 1,0^^11,-1; 1,1) = (A a -A h )/A. 

They describe the splitting of the hyperfine states that intersect in the absent of A^ nsec . 
In order to find out the picture of the anticrossing process we shall choose then as basis 
the next four states of the electron-nuclear spin system with M + m = —1 : \S, M; I, m) = 
|1, — 1; 0, 0), |1, — 1; 1, 0), |1, 0; 1, —1) and |0, 0; 1, —1). The reduced Hamiltonian can be rep- 
resented then by matrix 



( J/A-2fi B B {A a -A h )/A -(A a -A h )/A (A a + A h )/A \ 

(A a -A h )/A J/4-2/i B S (A a + A h )/4 —{A a - A h )/A 

-{A a -A h )/A (A a + A h )/4 g^B + J/A —(A a — A h )/A 

V (A a + A h )/A {A a -A h )/A -(A a -A h )/A g^B-3J/Aj 



(36) 



If for simplicity to assume here A a = A\ } = A, the matrix eigenvalues can be expressed 
by the simple equation: 

(J/4 - 2fi B B - E) ■ (gxfj.xB - J /A - E) - A 2 /Af - 

- [(J/4 - 2fi B B - E) ■ J/2} 2 = 0. (37) 

In this case Hamiltonian is symmetric, as in the absence of hyperfine interaction, with 
respect to interchanging donors a and b and eigenstates of Hamiltonian are divided into two 



In we have made mistakes for the values AE\(0, 0; 0, 0) and A£i(0, 0; 1, 0). 
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symmetric (I,S — 1) and two antisymmetric (I + S = 1) states, that are correspondent to 
the superpositions of the paired states |1, —1; 1, 0), |1, 0; 1, —1) and |1, —1; 0, 0), |0, 0; 1, — 1) 



with eigenvalues |L6 



E s ± = g N ^B/2 + J/A-fi B B±J(fi B B + g N ^ N B/2y + (A/2y 



El = g N fi N B/2-J/4-fi B B±^/(fi B B + g N ^B/2-J/2y + (A/2) 2 . (38) 

The splitting of the states |1, —1; 1, 0) and |1, —1; 0, 0) degenerated in the first approxi- 
mation is now 

E s -E & = 2irhv 3 OV^L f or j ^ 2 a B B, 

J 2fi B B - J 2fi B B ' P ' 

E s _-E* = 2-kHvj ~ A/2, for J = 2fi B B. (39) 

The energy 2irhi , j is the energy of transition of paired nuclear spins from the upper 
triplet state to the singlet state. It is known that such splitting can be described by spin 
Hamiltonian like / a blalb ; where J a b = 2irhvj is the indirect "exchange" integral for nuclear 
spins. For B = 2T and J/27fh = 30 GHz < 2[i B B /2nTi = 57 GHz, z/j was estimated at 
i>j = 75 kHz < i/A = 92.6 MHz ||. Note that the frequency z/j grows as we approach the 
crossing point of the unperturbed levels E (l, —1) and E Q (0, 0) and vanishes at J =3- 0. 

As it follows from (|38|), the electron-nuclear antisymmetric state 1 1, — 1; 0, 0) for J < 
2fi B B undergoes the transformation into anticrossed antisymmetric state |0, 0; 1, —1) for J > 
2fi B B. That is, the electron and nuclear subsystems exchange their states with the opposite 
transfer of nuclear spin states near the crossing point. The electron-nuclear symmetric state 
1 1, — 1; 1, 0) for J < 2fi B B does not anticross with any ground electronic states 1 0, 0) for 
J > 2fi B B and does not change nuclear spin states. Note that the splitting of symmetric 
and antisymmetric states levels for J 3> 2fi B B 

E s _ — E't = 2rthvj ~ J - 2fi B B (40) 

is more then ( |39| ) for J <ti 2fi B B. Under these conditions, the constant of indirect scalar 
interaction 7 a b between the nuclear spins of the donor atoms a and b also radically changes. 

The hyperfine splitting of triplet and singlet electronic levels near the point C is shown 
schematically in Fig. 3. 



6 Measurement of nuclear spin states 

Measurement of individual nuclear spin states is one of the important problems. We will 
discuss here this very briefly. 

The nuclear spin state is set for J <ti 2fi B B, where nuclear spins are handled with NMR 
methods by applying RF pulses of resonance frequency. It is suggested to measure this state 
in two stages 0. Suppose that after quantum computation at J <^ 2ji B B two electrons, 
both on its donor atom, are initially in the triplet state that is in the |1, —1) state, the whole 
electron-nuclear system is either in |1, — 1;1,0) or in 1 1 , — 1 ; 0, 0) state. By adiabatically 
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increasing of the exchange parameter to J > I^B in a process of the crossing point passage 
we lead to the electron-nuclear system transition from one antisymmetric state |1, — 1; 0, 0) 
to other |0, 0; 1, — 1) at the same total spin projection and allow transfer the information 
from nuclear to the electron spin subsystem. 

In addition, if the energy of bound for an electron at one neutral donor (it is usually 
small) is more then its energy of attraction to the neighboring ionized donor {D + — state), 
the electron will be found near the neutral donor (D~— state or helium-like atom) and charge 
transfer from one donor to the other will occur. This may be reach also by the corresponding 
change of the gate A electric potential. Therefore, a charge transfer from one donor to 
another takes place. It is supposed that this process can be detected with highly sensitive 
single- electron capacitive techniques. 

It is important to note that the electrical measurements time should be significant less 
then the time of electron spin-lattice relaxation, that at low temperatures exceeds 10 3 s. In 
the case of an electrical control and electrical measurements the states of qubits the noise 
in electrical circuits are also significant. They cause fluctuations of the nuclear resonant 
frequencies v& and vj and lead in result to decoherence and dissipation of quantum states 
and to quantum errors in the process of calculations. 

The detailed theoretical investigation of the electrical single-spin measurement method of 
electron and nuclear spin states by using single-electron transistors in simpler semiconductor 
structures with one double donor (in particular Te) was given in [[17]]. Another scheme for 
measurement of the state of single spin based on the single-electron turnstile and injection 
of spin polarized electrons from magnetic metal contacts is proposed in [0 



7 Bulk-ensemble solid-state quantum NMR computers 

The approach suggested in || can solve a number of problems typical for the considered 
individual- access computing. First, to initialize the starting state of the system, it will suffice 
to go to low temperatures T <C 2^B/k, at which each electronic spin of a donor atom is in 
the pure quantum ground state and thermal fluctuations are suppressed. Accordingly, the 
nuclear spins are also in the pure state. No special procedure to separate pseudopure states 
is required in this case. Second, the number of spins-qubits in this case is unlimited, which 
opens the door to creating an actual many-qubit quantum computer instead of demonstrating 
only quantum principles with the simplest spin systems on organic molecules. Third, there 
appears a possibility to detect, with single-electron techniques, electron transfer between 
adjacent A gates and thus to sense the state of individual nuclear spins. In ||, the possibility 
of creating hybrid quantum-classical systems was also noted. Here, quantum nuclear-spin 
devices are supplemented with conventional integrated circuits when many-qubit systems 
are difficult to obtain. 

However, there are certain difficulties in implementing a quantum computer on individual 
donor atoms that was suggested in |J. These are, first of all, small signal from the spin of 
an individual atom, the need for use of high sensitive single-electron measurements, then the 
need for precision regular donor arrangement with nanometer scale that matches the gate 
chain, and at last the need for use of low temperatures. 
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As an alternative, we will briefly discuss here another feasibility of a bulk-ensemble silicon 
quantum computer. In this case, unlike the structure suggested in ||, gates A and J form a 
chain of narrow (I a ~ 10 nm) and long (several micrometers) strips along which donor atoms 
L distant from each other are placed (Fig. 4). Thus, they form a quasi-one-dimensional 
regular structure. It can be considered also a random distribution of donors under gates. 

We shall consider here only two cases. 

Let L be so much larger than I that exchange spin interaction between electrons of donor 
atoms along the strip gates (y— axis) is negligibly small; i.e., J(L), kT <C J (I), 2fi B B. Then, 
such a system breaks down into independent chains of the donor atoms in the direction 
transverse to the gates (x— axis). In that case, the regularity of donor structure along the 
strip gates does not play any role. The gates form an ensemble of independent equivalent 
multiple-qubit computers — artificial "molecules" whose electronic spins are initially aligned 
with the field. Accordingly, all nuclear spins-qubits are also similarly oriented. An output 
signal in this system, as in liquids, will be proportional to the number of the "molecules" 
or donor atoms under a strip gate A. For example, for number of "molecules" more then 
thousand the length of the strips should be more then ten micron. 

One would expect that with the bulk-ensemble approach, where many "molecules" work 
simultaneously, electrical measurements would be greatly simplified. 

A more complicated situation appears when L < I. In this case, J(L) ^> 2/xb-B, J (I) 3> 
kT, and exchange interaction between localized electronic spins along the strip gates will 
be favorable for producing one-dimensional antiferromagnetically ordered chains, that may 
be considered at low temperatures as a pure macroscopic quantum state. The neighboring 
electronic spins in the antiferromagnetic chain in the absence of the field are oppositely 
oriented. Due to hyperfine interaction in the ground state, nuclear spins will also be oriented 
according to the electronic spin direction in the resultant field and will form in chains an 
antiferromagnetic type order. 

The nuclear resonant frequencies u\ of neighboring nuclear spins will be different for each 
of the magnetic one-dimensional subarrays in the chain: 2whu^ ~ Ig^fi^B ± 4|. As in the 
case of electronic spins oriented only in the external field, the RF field has a greater effect 
on a nuclear spin because of hyperfine interaction. A resulting NMR signal tuned only to 
one of the frequencies will be proportional to half the number of nuclear spins in the 
antiferromagnetic chain or half the number of donor atoms under the gate A. 

The selective RF it— impulse tuned to one of the resonant frequency can invert the nuclear 
spins state of one subarray in a chain and the nuclear spins of whole chain will acquire 
quasi-one dimensional ferromagnetic ordering while the neighboring chain will remain in the 
antiferromagnetic state. 

Let now the electronic spins of the such two neighboring chains are also setting for 
J (I) <C 2/ib-B, that is let they and half the number of nuclear spins in the chains are in the 
state with the same orientations (triplet states of the spin pairs), while other half the number 
of nuclear spins is in the state with opposite orientations (singlet states of spin pairs). By 
adiabatically increasing of the exchange parameter by means of J— gate potential varying to 
J (I) > 2/iB-B in the passage of the crossing point we lead for each electron-nuclear pair of two 
neighboring chains to the transition from one antisymmetric state into other antisymmetric 
state. The triplet electronic pairs ground state of two neighbor chains must pass to the singlet 
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state. The subarray of nuclear spins which is in the state with opposite orientation of nuclear 
spins relative to the subarray of neighboring chain will transit in the state with the same 
orientations. The electron subsystem of two chains simultaneously passes in antiphase state 
and at the same time to the half period shift of antiferromagnetically ordered electronic 
spins in the one chain relative to the other. The nuclear subsystem of the both chains 
becomes ferromagnetically ordered. Again, the Pauli principle allows two electrons to be 
now on the one site, if it is energetically profitable, that it allows transfer electrons between 
the neighboring chains and destruction the antiferromagnetic ordering in the chains. This 
process can be detected more easy than by means of highly sensitive single-electron spin 
states detecting methods. 

The relaxation time T 2 of nuclear spins, which characterizes the decoherence time of their 
quantum states, depends largely on fluctuations of local fields that are determined mainly by 
interactions between nuclear and electronic spins of different atoms. If the times of electronic 
spin-spin and spin-lattice relaxations in a semiconductor (r 2 and n, respectively) are short, 
that play the correlation times of fluctuating fields, electron-nuclear dipole-dipole and scalar 
interactions will be strongly averaged, like in liquids, and the time T 2 will be large. Only 
the scalar part of indirect interaction between the nuclear spins of donors J a bI a Ib will remain 
non-averaged. 

However, when the states of individual nuclear spins are determined by electrical tech- 
niques, the time t\ of electronic spin-lattice relaxation must be sufficiently large (several 
thousands of seconds) for the state of electronic spin have no time to relax during the elec- 
tric measurements. This again indicates the need for low temperatures || [l6j . In that case, 
for large enough distance between donors, large times T\ and T 2 are provided mainly by 
the small amplitudes of fluctuating local hyperfine fields. In insulators at low temperatures 
this fields are determined by electron-phonon mechanism. In particular it was demonstrated 



IS] that the phonon induced spin- lattice nuclear relaxation time T\ for low temperatures in 



silicon is very large. In the case of antiferromagnetically ordering at low temperatures the 
fluctuating local fields are determined also by the interaction of nuclear spin with electron 
spin waves 
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The decoherence is associated also with noise voltages across the gates A. To suppress 
the decoherence of different origin ideas similar to those underlying high-resolution NMR 



methods may be used pO 



Upon forming structures with closely spaced and narrow gates A, it would be more 
appropriate the lower arrangement of strip gates J. In this case, gates A can be made 
wider for the same interdonor distance along the x— axis. If we completely do away with 
electrical measurements of individual states, the control gates J become unneeded and the 
semiconductor structure simpler. Scalar interaction between nuclear spins of donor atoms 
through their electronic states, as well as their resonant frequency, will be controlled only 
by gates A and external magnetic field. 

An another interesting approach to implementation of a solid state ensemble quantum 
computer was described in [^2] , where it was considered a lattice of nuclear spins 1/2 having 
periodic structure of ABCABCABC. . . type in two or three dimensions, where A,B,C-are 
nuclear spins only of three types with different resonance frequencies. It is supposed that 
the nuclei at three sites are embedded in a crystal lattice of some solid state compound with 
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non-spinning nuclei and initialized with all spins at the ground |0) state. Each ABC-units 
of this superlattice can be used to store quantum information by setting one of spin up or 
down. This information can be moved around via some quantum cellular shifting mechanism. 
Cascading unitary quantum SWAP operations of A-^> B, B<^ C, C<S=^ A, A<^> B, . . . , every 
which is achieved by cascading three quantum controlled NOT gates, can be used for this 
process. An ancillary donor nucleus D with spin 1/2 in the proximity of an A-site can serve 
as the input /output port. 

Universal quantum logic is implemented with the aid of two-body interactions between 
two spins at D and a nearby site. A local environment region, about tens of ABC-units 
along all three dimensions provides a large quantum system with a number of qubits over 
thousands and only three types of nuclear spins. The whole crystal contains a huge ensemble 
of such identical NMR quantum computers — large artificial "molecules". The information 
in this case can be input by setting the D-spins to desired states then transferring to the 
nearest A-spins via the quantum SWAP operation. All these can be done selectively if the 
specified quantum transitions are driven by RF fields with distinguished frequencies. After 
this, the D-spins are reset to the |0) state. 

The state of any qubits can be measured by moving it to the A-site nearest to the D- 
spins, then by swapping A=>- D. Finally the state of the D-spins is measured by using NMR 
techniques. The proposed NMR quantum computer by the authors opinion can be easily 
scaled up and may work at low temperature to overcome the problem of exponential decay 
in signal-to-noise ratio. However, authors of the article |22| did not made any numerical 
estimations and discuss any possibility of a concrete realization of this idea. 

In conclusion, we will consider briefly possible nonelectric approaches that would allow 
ensemble solid state computers to actually operate at room temperature. Several such ap- 
proaches are available from the literature. 

Highly sensitive optical measurements of a change in Overhauser shift in exciton optical 
spectra were discussed in [|23|| . The shift changed as a result of nuclear polarization in 
GaAs/Al x Ga 1 _ x As quantum dots. In-plane spatial resolution of ~ 10 nm and a sensitivity of 
10 4 nuclear spins at 6 K were achieved. This ensemble variant, based on optical detection, 
requires considerable redesign of the structure for laser radiation can strike required regions. 
Note that the sensitivity attained in [^| is sufficient just for the ensemble approach. However, 
there are not well-known sensitive optical measurements for semiconductor structures in 
which silicon, germanium etc. are used. 

At higher temperatures, a number of spins will pass into the excited state, while others 
will remain in the ground state. Of the latter spins, one can construct the pseudopure state, 
using methods mentioned in the Introduction for the ensemble approach. Another way to 
keep as many nuclear spins as possible in the ground state is to orient them not only with a 
permanent external magnetic field but also with dynamic methods, such as optical pumping 



2J], which make possible the orientation of nuclear spins as high as nearly 100% even at 
room temperature. However, this can be done only for direct-gap semiconductors like III-V 
compounds. In the case of silicon or germanium this property of the band structure may be 
received, if the single crystals will be properly strained. 

The authors are grateful to V.A.Kokin and L.E.Fedichkin for critical reading of the article 
and for useful remarks. 
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Fig. 1. Two cells of the semiconducting structure with donor atoms a and b. /a 
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Fig. 2. Energy levels in the electronic-nuclear system of an individual donor atom. 
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Fig. 3. Energy levels of a dual-spin system subjected to a magnetic field in the absence 
of nuclear spins. The inset schematically shows hyperfme splitting of the ground electronic 
states near the crossing point C for A a ^ A\>. The circles indicate the anticrossing points of 
hyperfme states. 
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Fig. 4. Fragment of the suggested ensemble structure (region only under gate A). 



21 



